function y = Cauchy_PDF(x,m,b)
y = (1/pi) * (b ./ ((x - m).^2 + b.^2));
end